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Abstract 

The vector form factor f+{t) of the semileptonic decay D Kii', measured re- 
cently with a high accuracy, can be used to determine the strong couphng constant 
9d*dk- The latter is related to the normalised coupling g releveant in heavy- 
meson chiral perturbation theory. This determination relies on the estimation of 
the residue of the form factor at the D* pole and thus on an extrapolation of the 
form factor in the unphysical region {mo — t^k)'^ < t < {niD + mK)'^ ■ We test this 
extrapolation for several parametrisations of the form factors by determining the 
value of g, whose value can be compared to other (experimental and theoretical) 
estimates. Several unsophisticated parametrisations, differing by the amount of 
physical information that they embed, are shown to pass this test. An apparently 
more elaborated parametrisation of form factors, the so-called z-expansion, is at 
variance with the other models, and we point out some significant shortcomings of 
this parametrisation for the problem under consideration. 



The weak transitions from one meson to another provide very interesting tests of 
our understanding of the Standard Model both in the weak and strong sectors. On 
one hand, elements of the Cabibbo-Kobayashi-Maskawa matrix can be determined by 
comparing the experimentally measured decay rate with a theoretical value obtained at 
one or several kinematic points: for instance B — >■ Diu for \Vcb\, B — > irlv for \Vub\i or 
K nil/ for \Vus\- On the other hand, the shape of the spectrum provides a stringent 
check of our description of the dynamics of hadrons governed by QCD. 

The hadronic physics of the problem is encoded in form factors, which are complicated 
nonperturbative objects. Due to our limited theoretical knowledge of these objects, 
many parametrisations are available, trying to include as much information as possible 
on the (known or assumed) dynamics of the corresponding mesons. Recently, the BaBar 
collaboration has performed a very accurate analysis of the K~e^Ve decay [1] 

(see also the work done by the CLEO collaboration [2J). The relevant form factors are 
defined as 

rnvTsiAm) = {p,+p', - %'^^^^^) hit) + g.'^^^^^m (i) 

with q = p — p' and t = q^. Ref. [1] determined very precisely the t-dependence of the 
vector form factor /+ (the scalar form factor /q could not be obtained because of the very 
light mass of the electron). We can provide a general expression of this form factor in 
terms of an unsubtracted dispersion relation (such an unsubtracted dispersion relation 
is allowed by the asymptotic behaviour expected from QCD): 

/,(i) = 5|5(M + if,,5!LM:) (2) 



mi,, —t Tx Jt+ t' — t 

^ s 

where 

t+ = (mi5 + mxY t_ = [rriD - niKY (3) 

meaning that we have a pole at m\), (m/). = 2.112 GeV) and a cut from the DK 
continuum 0. We have defined Res(/+) so that it be positive. 

The physical region for the semileptonic decay is < t < t_, where high-precision 
experimental information is available pLj. Although it corresponds to an unphysical 
point between the semileptonic region and the DK cut, the pole residue is physically 
interesting. Indeed it is related to the D-K-D* strong coupling (with a subsequent 
weak decay of D* into a lepton pair), which itself hnked to go'Dn by SU{3) symmetry. 
This latter coupling is a fundamental ingredient for Heavy Meson Chiral Perturbation 
Theory and allows one to compute the effect of pion and kaon exchanges on the 

non-perturbative dynamics of heavy-light mesons (see refs. [7f8|9] for recent applications 
to semileptonic decays and other processes, as well as pending issues in this field). 



^In principle, the cut could begin at a much lower value, i.e. DgT^^, near the D*. Indeed, the 
DsTt'^ state has the appropriate quantum numbers to contribute to the form factor. However, DsT^^ — > 
DK scattering is exotic and there is no resonance able to yield a sizeable amplitude for this process. 
Therefore, this cut is likely tiny below the DK threshold, and we will not consider it in the following. 
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We want to exploit the good experimental knowledge on /+ in the physical region to 
determine the residue Res(/+) and thus gD*DK by an extrapolation of the form factor 
from i_ to m\,* : 

m = ^sisi2El^ + ... (4) 

where is the decay constant relative to the purely leptonic decay of the D*. gD*DK 
is related to another significant quantity, the normalised matrix element of the axial 
current between the B and B* mesons denoted g. 

For this extrapolation, we can devise analytical models or parametrisations describing 
the physical region as a function of t, and extrapolate the expression into the region near 
the pole. We then deduce the residue according to: 

Res(/+)- hm {ml,-t)U{t) (5) 
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It proves interesting to define a function Res(/+)(i), which we call the residue function, 
according to: 

Res(/+)(t) = {m% - t)U{t) Res(/+) = Res(/+)(m^,) (6) 

We can hope to get reasonable extrapolations by considering smooth representations of 
Res(/_|_)(t). Since we have factorised the pole denominator, the residue function has 
singularities only due to the cut. 

Models will show differences in the way that they represent these singularities. We 
propose to use the value of ^ as a test to select the most appropriate models among the 
available ones. Indeed, we have independent experimental and theoretetical information 
on this coupling. Models with reasonable physical assumptions (i.e. the singularities 
along the cut) should be able to provide reasonable extrapolations from the semileptonic 
region to the pole, and thus to yield values of g in good agreement with our current 
knowledge. 

Of course, the various heavy-to-light semileptonic decays could yield similar strong 
couplings, somewhat related by heavy quark or SU (3) symmetries. To evaluate the 
specific interest of the reaction D — > Klu in this respect, we can notice the following 
facts. Since 

i_/m^j = 0.42 t+/m]j. = 1.24 (7) 

the Dl pole is relatively close to the threshold of the cut for DK production, but far 
from the physical region for the semileptonic decay; in addition, the physical interval 
is not very large. This can be compared with other semileptonic heavy-to-light decays. 
B,D ^ Triu, exhibit the advantage of a pole located very close to the physical region, 
and specifically for B Trii/, the physical region is very extended. However the cut is 
now so close to the investigated pole that it may cause problems. For B — > ulu, the 
smallness of \Vub\ makes an accurate measurement of the form factor near the t = t_ 
endpoint very difficult: the decay rate is very small at this endpoint where its value is 
crucial for a good extrapolation. 
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D Tilv seems more promising in this respect, and it is complemented interestingly 
hy D ^ Kiu. In the latter case, the pole lies farther than for D tt/z/, but it is better 
separated from the cut, and it is easier to observe because of the large value of the CKM 
matrix angle Vcs (compared to the Cabibbo-suppressed D — > vr transition). On one 
hand this could be considered unfavourable: the extrapolation will remain undoubtedly 
uncertain, given the limited precision of the data, since many models may be very close 
in the physical region, but will greatly differ near the pole. On the other hand, it 
provides an interesting opportunity to test these parametrisations in a situation where 
their differences will be enhanced and thus easier to discuss. 

This note is organised as follows. In Sec. 1, we discuss the essential ingredient of the 
dispersive representation of the form factor, namely the DK cut. In Sec. 2, we introduce 
several parametrisations which differ mainly through their approximate description of 
the cut. In Sec. 3, we discuss the extrapolation of the form factor in the unphysical 
region according to these parametrisations. In Sec. 4, we collect the resulting values for 
the hadronic coupling constants Qd'DK and g and comment on the discrepancies induced 
by the different parametrisations, before drawing a few conclusions in Sec. 5. 

1 The DK cut 

Theoretical inputs concerning the cut associated with DK production are essential since 
it governs the behaviour of the residue function Res(/+). There are none compelling. 
Yet we can consider two partial and complementary contributions from the low-energy 
DK continuum and from resonances. 

First, consider the contribution to the cut stemming from the DK low-energy contin- 
uum. In the elastic regime, between the DK threshold and the first inelastic threshold, 
one has 

Im/+(t) = VXtf['^it)Uit)eit - (Mn + MkY) (8) 

where A = (t — — tJ) = Atq^ and f^^ is the P-wave amplitude for DK DK 
scattering. Close to the DK threshold (small positive g^), the form factor has the 
following behaviour in 

f^^{t) ~ A,q' (9) 

so that 

Im /^^(t) ~ A,q' (10) 

which means that the imaginary part departs slowly from zero above the DK threshold: 
the cut should be very smooth. 

Secondly, there are resonances all along the cut, corresponding to the successive 
radial excitations of the D*, the first of which should be close to the DK threshold. 
Indeed, with an excitation energy of order 0.5 GeV added to the mass 0, we get a 
first radial excitation around 2.6 GeV, while the cut starts at -|- rriK — 2.4 GeV. If 

^This order of magnitude for the excitation energy is found in most quark models, see for instance 
ref. [ini discussed below. 
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we worked in the large- A^^c approximation of resonances with vanishing widths, the cut of 
the form factor would be largely dominated by the pole of the first excitation. On the 
other hand, in the actual world, most of the resonances acquire a broad width through 
their strong decays, will overlap and interfere, yielding a rather smooth cut. This is 
presumably true even for the lowest radial excitation, if the coupling constant is not 
exceptionally small, because of the allowed phase space, so that the impact of this first 
resonance on the cut is presumably mild. These arguments are supported by the recent 
observation of the Dsj{2700)^ by the Belle collaboration [llj. This meson is interpreted 
as a radially excited cs state with = 1 with a mass M = 2708 ± 9t\l MeV, quite 
far away from the start of the cut, and a fairly broad width T = 108 ± 23^i^ MeV. 

Even though we have arguments supporting the smooth rise of the cut above thresh- 
old, it would be useful to estimate the absolute strength of the cut in order to constrain 
our extrapolation more tightly. Ideally, a determination of the vacuum-to-Di^ matrix 
element would provide the basis for an interpolation between the physical region for the 
semileptonic and the region of the cut, rather than an extrapolation. Since this piece of 
information is currently lacking, we have to rely on parametrisations of the form factor. 

2 Parametrisations of the form factor 

We speak of "parametrisations" because we have admittedly little theory behind them. 
They are functions of t with the following main merits: 1) in the physical region they 
collect the data of the form factors in a both simple and accurate description; 2) they 
incorporate in a rough way the few safe statements that we can formulate about the 
form factors outside the physical region, so that they may still be trusted at least at 
a certain distance from the physical region (of course, it remains an assumption that 
we can trust them as far as the pole). As summarised in eq. ([2]), we know that there 
is a pole at the known value t = m'jj,, and a cut begins at t = {mo* + ttikY. The 
various parametrisations propose explicitly or implicitly different simplifications in the 
treatment of the cut, and we will sketch the salient features of the most widely used 
ones. 

In the case at hand, D — ^ Kiv, the experimental data is so accurate that the best 
option consists in fitting the parameters of the parametrisations on the data in the phys- 
ical region. It turns out that most of them fit the data equally well in the semileptonic 
region. A particular representation would exhibit a decisive advantage if it were able 
to describe the unphysical region between the semileptonic region and the cut with the 
same parameters. 

2.1 ^-expansion 

A parametrisation which is by now rather popular, yet apparently sophisticated, is the 
z-expansion [T^fT^] . which dates back to Boyd et al. (see ref. [H] and the works with 
Grinstein and Lebed quoted therein) where it was applied to heavy-to-heavy transitions. 
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In principle, it has been devised to offer a good representation of form factors and we 
recall some steps of its derivation below. 

One aims at writing down /+ as a series S aiz{ty . The expansion parameter z{t) is 
defined as: 

with a free parameter to- z{t) remains small over the whole physical region, and maps the 
complex t-plane into a disk of radius 1, with the cut in t being transformed into the circle 

\z\ 



1. Ref. [12] advocates the choice toA+ = 1 ~ y 1 ~ {to/rnjj, = 0.23), which 

implies that z{t) = occurs for t in the middle of the physical region [z{0) = 0.0515 and 
z{t-) = —0.0515]. The change of variables from t to z has often been used in the hope 
of getting a quicker convergence of the resulting series in the physical region. 

It has been proposed to combine this idea with analyticity constraints in order to 
constrain further the coefficients of the 2;-expansion [T2lll3fl4j . The starting point is the 
correlator 

U>"'{q) =ijd^x e*^"(0|T(s7'^c)(x) (s7''c)t(0)|0) = (g^g" - g^'')U^ (q^) + g^^U^iq^) 

(12) 

which defines two polarisation functions (transverse and longitudinal). First, these func- 
tions can be bounded approximately at large momenta using the Operator Product 
Expansion. Second, these functions can be expressed in terms of their imaginary part 
through a dispersion relation. Using unitarity, one can express this imaginary part as a 
sum of various (positive) contributions, corresponding to the squared modulus of matrix 
elements (iJ|s7^c|0), where H denotes any arbitrary single- or multiple-particle state 
with the appropriate quantum numbers (each contribution contains a a factor coming 
from the phase space). It means in particular that the imaginary part is larger than the 
sole contribution from H = DK, i.e. from /+, multiplied by a phase-space factor. 

By choosing multiplying the series in z by a carefully chosen factor denoted $, one can 
convert the bound induced on /+ by unitarity and OPE into a bound on the coefficients 
of the series in z. The ^-representation can be written as follows: 

, /+(t) _ P(0)$(0,to) l + f^.(t,to) + g.^(t,to) 

P-/+(0) Pm{t,to) l + f^z{0,to) + f^z^{0,t,) ^ > 

where 



$(t, to) = N{t+-t){^U - t + ^)-'i^U - t + ^U - to)i'sJt+ - t+^t+ - t^f^ (14) 

where is a numerical normalisation factor. It is unimportant here, since we consider 
only ratios of the form factor to its value at q^ = 0. The factor P(t) = 2;(t,m|,,) 
includes the D* pole. As explained above, the coefficients a„ are bounded by unitarity: 
J2 I flip < 1 with a properly chosen normalisation constant N p^l3lll4] . 

In heavy-to-light processes, this expansion seems less useful, because it appears that 
all the coefficients that can be determined are anyway much smaller than the unitarity 
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bound in absolute value [T], and nothing safe can be said about the relative magnitude 
of the coefficients with respect to the first one (see the fits in sec. 3). There remains a 
useful bound on the rest of the series, obtained by using Cauchy's inequality: 
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2 



(15) 



This bound, given the observed values of the form factor, is sufficiently small to be 
meaningful: it shows that the necessary number of terms should be about 2 or 3 to get 
a one percent accuracy in the semileptonic region. Actually, it was found that 2 or 3 
terms are sufficient to fit the data pQ, the last coefficient (02) being already affected by a 
large uncertainty. Unfortunately, this accurate representation in the semileptonic region 
does not warrant an accurate extrapolation. Indeed, in the region t_ < t < t+, z{t) by 
no means remains small in this region t^ < t < t+, since at the pole, z{mj^*) = —0.340 
and at the beginning of the cut, z(t+) = —1. In this region, the bound eq. (ITSl) becomes 
ineffective, or indicates that very many terms could be necessary. 



2.2 Alternative representations 

This parametrisation can be compared with other less sophisticated parametrisations 



Pole : 



BK 



Linear : 



Quadratic : 
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(16) 
(17) 

(18) 
(19) 



The pole parametrisation is certainly too naive, since there is no reason for the 
lowest lying pole to saturate the form factor. Not surprisingly, one finds the pole 
much below its actual position (the D* mass). We will not use this parametrisation 
to perform an extrapolation of /+ to the D* mass; it serves only to show that the 
achieved experimental accuracy requires to go beyond the common assumption of 
the dominance by the lowest pole in the physical region. 



The BK parametrisation is well known in the analysis of -B-decays, and it is 
mainly motived by the scaling laws of the form factors in the heavy quark limit. 
For a D decay, this motivation is certainly weaker, but it is useful by providing a 
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crude representation of the cut by an additional pole, with a free mass which comes 
out neatly above the start of the cut, at t = m£,*/a, with a smaller than 1 to be 
consistent with the location of the cut. Let us stress that this second pole is not to 
be identified with any particular resonance. This effective pole sums up the effect 
of infinitely many resonances with positive or negative contribution (note that the 
additional BK pole has a contribution opposite in sign to D*). 

• Almost equivalent in practice, the "linear" and "quadratic" parametrisations follow 
old ideas from current algebra. They assume that the amplitude have a polynomial 
dependence of low degree in the momenta, once the lowest lying (ground state) 
poles have been factored out. Such an expansion is justified in the case of a weak 
cut, which can be approximated by its expansion in t with a reasonable accuracy. 

• Another parametrisation can be obtained by considering the z-expansion and re- 
placing $ with a constant. Obviously, no unitarity bounds can be derived in such a 
case: we have just performed a change of variable (reexpressing a t-series into a z- 
series) and extracted the pole. In order to distinguish these two versions of the 
2;-expansion method, we call "simple ^-expansion" the parametrisation with $(t) 
set to 1, and "unitary 2;-expansion" the expansion relying on eq. f|T^ . which allows 
to exploit unitarity constraints (at least in principle) and described for instance in 
refs. IT21IT31IT1. 



3 Fit of the parametrisations to D ^ Kiv data 

Before extrapolating these various parametrisation to the pole, we need to determine 
their parameters from the data in the physical region. A fit of z-expansion, BK and 
pole parametrisations has been made in ref. pQ with data corrected for radiative effects. 
Unfortunately, the data (values of the form factors and correlation matrix) provided 
in this reference are only given before the correction of radiative effects. According to 
ref. the main correction corresponds to an increase of the value of the first bin. 
Therefore, for each parametrisation we can compare three different fits depending on the 
set of data: 

• The data before radiative corrections and the correlation matrix before radiative 
corrections (performed by us), 

• The data with the first bin increased and the correlation matrix before radiative 
corrections (performed by us), 

• The data and the correlation matrix with full radiative corrections (performed in 
ref. [1] and by us). 

The parameters of the different models are collected in Table [H Apart from the pole 
model, which does not fit the data very well, the other models yield very similar fits. 
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Mo del 


Model param 


No radiative corr 


First bin corrected 


Full radiative corr 


Unitary z-exp 


ai/ao 




— Z.4o ± U.Zi3 


-1- n ofi 
— z.o it U.Zo 




02/00 


— 4. ZD it 0.0 ( 


— z.yo it 0.0 ( 


U.o it 1.0 




correl 


-U.oz 


-U.oz 


-U.ob 


JDIV 


a 


U.oo ± U.Uo 


U.o / it U.U4 


n Q77 -I- n nQ7 
U.o M it U.Uo 1 


role 




1 _i_ n no 
i.oy ± U.Uz 


1 on _i_ n no 
i.oy ± U.UZ 


i.oo4 it U.Uiy 


Linear 


Cl 


0.43 ± 0.04 


0.43 ± 0.05 


0.42 ± 0.05 


Quadratic 


Cl 


0.45 ±0.15 


0.41 ±0.15 


0.33 ±0.15 




C2 


-0.036 ± 0.42 


0.046 ±0.42 


0.25 ±0.42 




correl 


-0.96 


-0.96 


-0.96 


Simple z-exp 


ai/ao 


-4.27 ±0.26 


-4.28 ±0.26 


-4.34 ±0.26 




02/00 


3.13 ±6.73 


4.41 ± 6.74 


7.78 ± 6.79 




correl 


-0.85 


-0.85 


-0.85 



Table 1: Parameters of the different parametrisations of /+ from a fit to D — ^ Kiu data, 
with various treatments of radiative corrections (not included, partially included, fully 
included) . 

The corresponding minimal values of the remains between 7 and 8, for 7 or 8 degrees 
of freedom (fit with 2 or 1 parameter). 

The first three entries of the last column are taken from ref . [TIITB] , where we combined 
systematic and statistic uncertainties in quadrature. Unfortunately, the three other 
parametrisations (linear, quadratic and simple ^-expansion) have not been considered 
in [T|. As a poor man's way of getting some information on the impact of radiative 
corrections for these two parametrisations, we have performed a fit with: 

• The correlation matrix provided in ref. [1] 

• The central values in each bin obtained from the 2;-parametrisation, using the 
values of ai/ao and 02/00 obtained including the radiative corrections. 

We have checked that this procedure leads to values and uncertainties for the fitted 
parameters of the different models which are very similar to those quoted in ref. [1]. 
However, this procedure is admittedly a very imperfect attempt of getting a handle on 
radiative corrections, and should be replaced by a full treatment of radiative corrections 
for these three parametrisations, which can be done only by our experimental colleagues. 

Figs. [Hand [2] illustrate the range of variation allowed for the residue function accord- 
ing to the different parametrisations, respectively before and after radiative corrections 
are taken into account. A striking difference can be observed between the (unitary) 
z representation and the other parametrisations, the errors in the first case exploding 
beyond t/m?j^, ~ 0.8. This behaviour can be understood from the explicit expression of 
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Figure 1: Range of variation for the function Res = f+(t){m'jj* — t) related to 

the residue of the form factor at the D* pole as a function of t/Ml,,, for the different 
parametrisations considered. No radiative corrections are included. 



the z-expansion. One must note the important factor t^. — t in the definition of $(t,to), 
eq. ( |T^ . This factor implies that the z-expansion includes a spurious, unwanted pole at 
the threshold of the cut, a feature which is very disadvantageous when the expansion is 
used not far from the cut, as needed to make our extrapolation. To be used there, the 
z-expansion requires in principle a very large number of terms in the series S a„2;", and of 
course so many a's cannot be determined from the knowledge of the physical region. At 
most, the first two or three coefficients can be estimated, so that the spurious threshold 
pole has a very strong effect on extrapolation. At the start of the cut, the 2;-expansion 
becomes meaningless. 

The complicated form of $ was chosen to translate the unitarity conditions on /+ 
into a constraint on the coefficients a„ of the 2;-expansion. In particular, the troublesome 
threshold pole is related to the two-body phase space arising in the unitarity bound on 
/+. On the other hand, it turns out that the bound on the coefficients a„ is far from 
saturated, and thus of little impact on the problem at hand, so that it can be left out 
without consequence, as exemplified by the "simple" 2;-expansion (with $ = 1). 

It is also interesting to compare this parametrisation with the BK parametrisation. 
Both representations have a pole in addition to D*, but the former fixes it at threshold 
while the latter keeps this pole at t = m^i./a as a free parameter. This effective pole 
sums up the contributions of a tower of resonances to the cut. We find that a < 1, 
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Figure 2: Range of variation of the function Res f+(t) = f+{t){m1), — t) related to 
the residue of the form factor at the D* pole, as a function of t/M'jy^,, for the different 
parametrisations considered. The effect of radiative corrections is included according to 
ref. |T] in the case of unitary 2;-expansion and BK models (see the text for the treatment 
of radiative corrections for the other models). 



consistently with the location of the cut. Let us also recall that one should have a ^ 1 
in the heavy-quark limit and large-energy release limit: therefore, it is consistent to get 
something rather different from 1 in the decay D Ktv ^ 



4 The strong coupling gD*DK 

Extrapolating the results of the previous section to the pole and propagating the 
errors, we obtain Qd'DK through eq. (jl]), which is related to the matrix element of the 
axial current between the B and B* mesons. After a suitable normalisation of states, 

•^The BK type of parametrisation can be used, and has indeed been used prior to the recent ex- 
perimental measurement, in order to represent model or lattice QCD calculations of 13 ^ Klv in the 
physical region, and deduce the residue at the pole. The authors of refs. [17ll8j found a — 0.24 from a 
quark model, and a lattice QCD computation [TO] led to a = 0.27(11) or 0.43(12), to be compared with 
the present 0.38. 
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this matrix element is denoted as g and one has 



9 = 




(20) 



g is finite in the heavy-quark mass hmit, and thus should vary slowly with the heavy- 
quark mass. Another important property is that the strong couplings of the K meson 
(Goldstone boson) with the D and any other Dg-type state satisfy an Adler-Weisberger 
sum rule, with all the intermediate states contributing positively, yielding a bound on 
the go'DK coupling. In terms of g, this bound is quite simple: 



The upper bound cannot be reached because there are indeed transitions from the D me- 
son to other states than the D*. The limit 1 corresponds to a completely non-relativistic 
calculation (for more details on g, see for instance ref. [lOj). 

Numerically, we take fx = 154 MeV, = 1.864 GeV, m^. = 2.112 GeV. Since 
fo* has not been measured directly, we must rely on lattice calculations. Np = 
calculations from the Rome [SU] and UKQCD [2T] lattice collaborations are available, 
with sizably different numbers. Unfortunately, neither papers quote their value of Jd* 
explicitly. In the case of the Rome collaboration, we have to pass through the ratio 
/^,//^. = 1.10 and the value of fo- = 258(14) (6) MeV. In the case of UKQCD, we have 
to convert their result from their dimensionless definition of fo* = 8.3 to our definition 
I'd* = "^D*/ I'd*^^'^^- In addition, it must be noticed that the mass mo* must be taken 
consistently from the lattice data at the same f3 and with the same way of fixing the 
lattice unit as fn*- We choose /? = 6.2 and we use fj, to fix the lattice unit (as done to 
compute fo*) so that rriD* = 2.061 GeV. Following this procedure, we obtain finally the 
central values /^j = 284 MeV (Rome) and /^j ~ 248 MeV (UKQCD). 

A very naive average yields the value fn* = 270 MeV that we use in the following. 
We do not quote errors on the auxiliary quantity fo* '■ it would represent a very involved 
task for a quantity that we use only as a reference scale to compare different parametri- 
sations. However, one should keep in mind that the value of this decay constant may 
be underestimated in view of the situation for the closely related quantity fo^- Indeed 
the same lattice groups predicted the latter around 230 MeV, while the most recent 
experimental average [22] yields 275 ± 10 MeV |. If fot were to be enhanced, gD*DK 
and g should be accordingly rescaled (and lowered). However, the hierarchy observed 
in the results of the extrapolations will not be affected by this change in the overall 
normalisation. 

We have discarded the simple pole model in the discussion of the residue: it does not 
account for the data properly, even with a flexible vector meson mass, and it would be 
meaningless to discuss the residue with a flctitious pole mass. 

* This relation is based on the heavy-quark expansion, but also on PCAC for the kaon and thus 
could be affected by corrections of order 30%. 

^The PDG Review of particle properties 2006 was quoting a still higher value: Jd^ — 294 ± 27 MeV. 



^1 < 1 



(21) 
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Model No radiative corr First bin corrected Full radiative corr 



Unitary z-exp 


16.2 ±7.8 


17.9 ±7.8 


23.0 ±8.7 


BK 


18.3 ± 1.1 


18.2 ± 1.0 


18.2 ± 1.2 


Linear 


16.3 ±0.6 


16.2 ±0.6 


16.1 ±0.6 


Quadratic 


16.1 ±6.5 


16.6 ±6.4 


18.1 ±6.4 


Simple z-exp 


15.4 ±3.2 


16.1 ±3.2 


18.2 ±3.1 



Table 2: Values of Qdidk obtained from different parametrisations of /+, with various 
treatments of radiative corrections (not included, partially included, fully included). 



Model 


No radiative corr 


First bin corrected 


Full radiative corr 


Unitary 2;-exp 


0.63 ±0.31 


0.70 ±0.30 


0.89 ±0.34 


BK 


0.71 ±0.04 


0.70 ± 0.04 


0.71 ±0.05 


Linear 


0.64 ±0.03 


0.63 ±0.03 


0.62 ±0.03 


Quadratic 


0.62 ±0.25 


0.64 ±0.25 


0.70 ±0.25 


Simple 2;-exp 


0.59 ±0.12 


0.62 ±0.13 


0.71 ±0.13 



Table 3: Values of g obtained from different parametrisations of /+, with various treat- 
ments of radiative corrections (not included, partially included, fully included). 

The results for Qdidk and g are collected in Tables [2] and [31 Models describing the 
singularities of the cut in an appropriate way should yield values of g in good agreement 
with what has been obtained from experimental measurement of strong decays or from 
the lattice calculations. 

The only indication from experiment is indirect because it concerns the SU (3) coun- 
terpart of gD*DK, 9d*Dit- Indeed, CLEO [23p^ has measured the decay width of 
D* — s> Z^vr, from which they estimate: 

go'D^ = 17.9 ± 0.3 ± 1.9 (22) 

leading to: 

^ = 0.59 ±0.07- 0.6 (23) 

The latter value should apply roughly for the D* — D transition through SU (3) flavour 
symmetry (the Dirac quark model and the lattice results suggest an increase when one 
increases the light quark mass up to the strange mass). 

Model approaches can be used to determine this coupling. The Dirac model, see for 
example ref. [ID], led to a result close to ^ ~ 0.6 for a static c quark, a number close to 
what was found later in the above experiment]^. The application of a dispersive approach 
to a constituant quark model led to g from 0.4 to 0.5 j28|IT8] . 

careful study based on QCD sum rules and including radiative corrections [3S] yields a much 
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The result of lattice QCD simulations for the D* — D transition is, by extrapolation 
to the u, d mass (ref. [29]): 



which is compatible with the experimental number quoted above. The lattice can in fact 
measure directly another SU{3) partner of the D* — D axial transition matrix element 
(corresponding to the Ai form factor, up to a mass factor very close to 1), i.e. D* — Ds, 
since the lattice direct measurement is close to the strange quark mass; the result is close 
to 0.7 from same reference, at f3 = 6.2. 

Comparing these values (all in the same range) with the above tables, which col- 
lect the coupling obtained by extrapolation of D ^ Kiu, we can make the following 
comments: 

• The result from the (unitary) 2;-expansion is sensitive to the extraction of the 
radiative corrections on the whole range of momenta. With a full treatment of 
radiative corrections, one obtains a neatly larger number than with the other 
parametrisations, and the central value of g is close to the upper bound set by 
the Adler-Weisberger sum rule. We interpret this as a consequence of the spurious 
pole at t = t+: this representation must be discarded when the unphysical region 
is concerned because of its unphysical behaviour at the threshold of the cut. 

• The four other results are rather close to each other which is encouraging, since 
they rely on varied but reasonable assumptions on the smoothness of the cut. In 
particular, the simple 2;-expansion yields a value in good agreement with the other 
models, confirming that the spurious pole in $, imposed by unitarity constraints, 
is the actual source of difficulties for the unitary 2;-expansion. 

• For these four models, the central value is in agreement with the independent 
determinations of g, in particular the experimental one. This provides further 
support for our preferred assumptions and parametrisations. 

• One can check a posteriori that the effect of the cut is varying very slowly with 
t. Let us take the linear parametrisation as an illustration. After subtracting the 
pole which is at its right position and strength, the remaining contribution to the 
form factor is just a constant, that is, the cut is represented by a constant c, with 
ci = — 0.42/+(0) (the negative sign indicates that it has a lowering effect on the 
form factor). With such a parametrisation, we are sensitive to the structure of the 
cut neither in the physical region - which is perhaps not too surprising - nor in a 
large part of the unphysical region t_ < t < t^. 

• Comparing the "quadratic" and "linear" fits, we observe that a two-parameter 
fit, which seems equally reasonable as a one-parameter parametrisation, yields no 

lower number. A possible way of curing this surprisingly small result has been proposed in [26j : on the 
hadronic side, a large contribution from radial excitations should be added to the standard perturbative 
continuum. Alternative estimations from light-cone sum rules for semileptonic B decays [27] point 
towards values of g closer to the ones collected here. 




(24) 
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change in the central value of g; however, the errors are increased leading to a 
loss of predictive power. This means that the parametrisations are not constrained 
enough by the data in the physical region for a compelling extrapolation. 

5 Conclusions 

We have exploited recent high-precision data of the BaBar collaboration on the D Kiu 
decay in order to test various parametrisations of heavy-to-light form factors. Indeed, the 
current theoretical description of such form factors remains largely incomplete, and it can 
be improved by a direct comparison with data. The accuracy of the experimental results 
yields good fits of these parametrisations of the vector form factor /_|_ in the physical 
region. But the value of the vector form factor outside this region is also of interest, since 
its residue at the mass of the D* meson is related to the strong coupling constant g. To 
extract this quantity, we need to extrapolate the form factors outside the semileptonic 
region, and thus to rely on the properties and assumptions of the various parametrisations 
for /+. On the other hand, the determination of hatg and its comparison with values 
from other approaches (measurements of strong decays, lattice computations) provides an 
interesting test of the parametrisations of the form factors, since it probes the differences 
in the way these parametrisations representation the physical singularities along the cut. 

In spite of the modesty of our approach, which is very phenomenological, some useful 
conclusions can be drawn: 

• The method of extrapolation using smooth models for the cut - either a remote 
effective pole (BK) or a low-order polynomial in t - has an encouraging success, 
since the D* residue appears quite compatible with values expected from very 
different considerations. The z expansion (in its unitary version, widely popular 
by now) is disfavored: it contains a spurious pole at the threshold of the DK cut 
which makes the extrapolation blow out of control in the unphysical region. 

• The uncertainty on fo* is still large; reducing it would in turn reduce that of the 
strong coupling constant and help to strengthen the previous conclusion. 

• The extrapolation has by itself a very large uncertainty, which cannot be reduced 
further if we know only the physical region. Some quantitative theoretical knowl- 
edge about the form factor on the cut , in particular at threshold (scattering 
length. . . ), would be a great help by transforming the extrapolation into an in- 
terpolation. Borrowing ideas from ref. [30], one could then use a sufficiently sub- 
tracted Omnes-Muskhelishvili representation to combine our knowledge on the cut 
near DK threshold and the physical semileptonic region. 

• Another improvement would consist in enlarging the interval of the fit. Such in- 
formation could be provided by lattice simulations computing the form factor for 
t < (indeed they can calculate the D — K matrix element without referring to the 
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D Klv transition). However, the improvement would be much milder, since the 
sensitivity to the cut decreases when one gets deeper into the region of negative t. 

The D Klv process is a particular decay among many similar semileptonic pro- 
cesses. It presents the advantage of a highly accurate knowledge of the form factor, 
and of a clear analytical structure, with the physical threshold, the pole and the cut 
neatly separated, which will not be the case for other processes. Our simple analysis 
of this particular process shows some interesting features of the different parametrisa- 
tions of the form factors currently used to analyse weak transitions from one meson to 
another. In particular, it should invite practitioners to proceed with care when they 
have to rely heavily on parametrisations to extract quantities of physical interest (such 
as CKM matrix elements or strong hadronic couplings) with a high accuracy. 
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